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BOREL THEOREMS FOR RANDOM MATRICES FROM THE 
CLASSICAL COMPACT SYMMETRIC SPACES^ 

By Benoit Collins and Michael Stolz 
Universite Claude Bernard Lyon 1 and Ruhr-Universitdt Bochum 

We study random vectors of the form (Tr(A(^'l/), . . . , Tr(A(''V)), 
where 1^ is a uniformly distributed element of a matrix version of a 
classical compact symmetric space, and the A'^"'^ are deterministic 
parameter matrices. We show that for increasing matrix sizes these 
random vectors converge to a joint Gaussian limit, and compute its 
covariances. This generalizes previous work of Diaconis et al. for Haar 
distributed matrices from the classical compact groups. The proof 
uses integration formulas, due to Collins and Sniady, for polynomial 
functions on the classical compact groups. 

1. Introduction. One way to extend the scope of random matrix theory 
is to revisit its roots in physics. Classically, random matrices were considered 
by physicists as toy models for operators associated to complicated quantum 
systems. The matrix ensemble was supposed to be compatible with the sym- 
metries of the quantum system. Analyzing the action of symmetry groups, 
Freeman Dyson [10] argued in 1962 that the good objects to study were 
the Gaussian Orthogonal, Unitary and Symplectic Ensembles (GOE, GUE, 
GSE) as models for Hamiltonians, and the Circular Orthogonal, Unitary 
and Symplectic Ensembles (COE, CUE, CSE) as models for time evolu- 
tion operators. The geometric fact behind Dyson's result is that the matrix 
ensembles in question are associated to classical symmetric spaces. 

For concreteness, let us review the case of GOE/COE: GOE(n) is a prob- 
ability measure on the space of real symmetric n x n matrices such that the 
matrix entries on and above the diagonal are independent centered Gaus- 
sian random variables. What is more, it is assumed that the variance of the 
diagonal entries is twice the variance of the off-diagonal entries. This guar- 
antees that the Lebesgue density of GOE(n) factorizes with respect to the 
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trace and is hence invariant under conjugation by orthogonal matrices. Let 
us give an abstract account of these weh-known facts: Write 

u := un := {X G C"^" :X' = -X}, 

the real vector space of skew-Hermitian complex nx n matrices, u admits a 
real linear involution 0{X) = X. Then we have the eigenspace decomposition 

u = t®p, 

where 

is the +1 eigenspace and 

is the —1 eigenspace. Note that 

ip = {XG]R"^":X' = X} 

is the space on which GOE(n) is supported. Applying the matrix exponential 
to the elements of u, we have 

hence e'^ is unitary. In fact, any unitary matrix arises in this way, and Un 
is the Lie algebra of the unitary group U„. t exponentiates to the special 
orthogonal group S0„ and is its Lie algebra. So, on this abstract level, 
GOE(n) is a probability distribution on ip, which is invariant under the 
action of the Lie group corresponding to 

For X G p, the image is both unitary and symmetric, and it turns out 
that 

exp{p) = {gg' ■.g£\Jn}=:S. 

For g G U„ write ^{g) = gg' ■ Then $~^(1) = O^- We can thus view $ as a 
bijection of U„/On onto S. The uniform distribution on Vn/On is Dyson's 
circular ensemble COE. Note that 0„ C is the set of fixed points of 
the involution ^(^f) = {g~^)' of U„. Coset spaces of this kind arise when 
Riemannian symmetric spaces (RSS) are represented as homogeneous spaces 
of Lie groups, and conversely, modulo some technicalities, coset spaces of this 
kind are RSS (see [14] for details). In this sense we can view S and p (hence 
ip) as representing an RSS. 

As a matter of fact, the matrix ensembles from Dyson's list are of Cartan 
class A, AI, All; that is, GUE(n) lives on iu, which is nothing else than the 
space of hermitian nxn matrices, GOE(n) is as above, and GSE(n) lives on 
the ip part in U2n, this time with respect to the involution 0{X) = JnXJ~^, 
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where J„ is as in (1) below. In recent years, it has emerged that the frame- 
work that led to Dyson's list is not flexible enough to meet the needs of 
condensed matter physics; see, for example, [1]. A classification program, 
in the spirit of Dyson, has been carried out by Heinzner, Huckleberry and 
Zirnbauer in [13], to the effect that all spaces ig, ip arise naturally as spaces 
of (toy models for) Hamiltonian operators in mesoscopic physics, where q is 
any classical compact Lie algebra and p is the —1 eigenspace of a classical 
"Cartan involution" of g (see Section 2 below). By the same token, all clas- 
sical compact Lie groups and spaces S are needed to model time evolution 
operators. 

This paper is devoted to the S side of the story. For the ig, ip side see 
[11]; for other aspects of the S side see [9]. We endow S with the uniform 
distribution and investigate the behavior of the matrix entries of a random 
element of S as matrix size tends to infinity. We obtain Gaussian limits 
for wide classes of vectors of linear combinations of matrix entries. The 
archetypical result in this direction is due to Emile Borel, who proved a 
century ago that for a uniformly distributed point X = {Xi, . . . , X„) on the 
unit sphere S"~^ C M", the scaled first coordinate -^/nXi converges weakly to 
the standard normal distribution as n— > oo. Further-reaching results about 
the entries of matrices from the classical groups have been obtained in [5, 6] 
and [15], the first reference containing illuminating remarks on the history of 
the topic. To prove the weak convergence results of this paper, we invoke the 
method of moments. To make it operational, we use what we term "extended 
Wick calculus." The key ideas can be traced back to Weingarten [17]. They 
have been developed by Collins [3] and Collins and Sniady [4]. There is a 
certain family resemblance between this calculus and the Diaconis-Evans- 
Mallows-Ram-Shahshahani approach to random matrix eigenvalues [7, 8], 
in that both exploit in different ways the invariant theory of the classical 
groups (see [16] and Section 3 below). 

The paper is organized as follows: In Section 2, we explain the basic 
setup and notation and state the main results. Then, in Section 3, we re- 
view extended Wick calculus and explain its connection with the Diaconis- 
Shahshahani integral. Finally, the main results are proven in Section 4. 

2. Notation and statement of the main results. For n G N let G = Gn 

be one of the unitary group Un, the (compact) orthogonal group 0„ and 
the (compact) symplectic group Sp2„. g' denotes the transpose of a matrix 
5, and g* denotes conjugate transpose. For n G N, write I„ € C'"^" for the 
identity matrix. For p, g, n £ N,p + g = n, set 



(1) 
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and 




Note that Ip^g = 1^^ and J' = J ^. We consider pairs {Gn,Oc), where 9c is 
an involutive automorphism of G„, given as follows: 



• For G, 



\Jsn, we set 



Oai{9) 
0Au{g) 




• For G, 



Osn, we set 



6'bdi(5') 
6'diii(5') 




• For Gn = Sp2„ , we set 

dcnig) = 'Ip,qg'Ip,q- 

Here, J = 5c = 2 for C = All, Dili and 5 = 1 otherwise. Write K = Kn = K^ 
for the fixed subgroups of the automorphisms Oq- Concretely, K^^ = 0^,, 



i^AII = Sp2„, K^''' = Up X U„ i^^DI ^ Op X 0„ i^DI" = C 02n, K^' = 
Un C Sp2„, if^" = Sp2p X Sp2g. 



Then the list of homogeneous spaces Gn/K^ is cum grano salis the list of 
classical compact symmetric spaces, as classified by Cartan. But it should be 
emphasized that we follow standard practice in probability theory in that we 
consider the full groups U„ and 0„ rather than SU„ and S0„, which would 
be the natural objects to consider in Lie theory. We will refer to {Gn,0c) 
as symmetry class C, and we will say that is of class A, is of class 
B/D and Sp„ is of class C. The differentials of the Oq are the involutive 
automorphisms of the Lie algebra g„ of G„ which were discussed in the 
Introduction. Write 



We will refer to $c as Cartan embedding and to Sc as the Cartan-embedded 
version of the compact symmetric space Gn/K^ of class C. To make the con- 
nection with the preliminary discussion of GOE/COE in the Introduction, 
note that the tangent space in 1 to 5c is the (—1) eigenspace p of the differ- 
ential of 9c- Adopting physics terminology, we call A, AI, All the "Wigner- 
Dyson classes," AIII, BDI and CII the "chiral classes" and all others "su- 
perconductor" or "Bogolioubov-de Gennes (BdG)" classes. Observe that the 



Mg) = ^Fi9) = gmg)) 



and Sc = S, 



C 



MGn) C G, 
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chiral classes are special in that they depend not only on the overall matrix 
size n (resp. 2n), but also on parameters p = p{n), q = q{n) with p + q = n. 

For completeness, define ^j^y'j-,, ^>[?^ as identity on \Jn,0n,Sp2n^ 
spectively. Write Xg = Haar measure on normalized as a prob- 

(n) 

ability measure. Write V = Vc = V^ for a random element of 5c, whose 
distribution is the push- forward of Xq under It is convenient to call 
classes A, AI, All, AIII the "A classes," B/D, BDI, Dili the "B/D classes," 
and C, CI, CII the "C classes." Set IK = C for the A-classes, IK = M for 
the B/D-classes and IK = EI for the C classes, EI denoting the skew field of 
quaternions. We embed the space EI"^" of quaternionic matrices into C^"^^" 
as 



X_ Y_ 

-Y X 



:X,y G 



-<nxn 



Recall that Sp2„ = M"^" n Ugn and 0„ = M"^" n U„. We endow IK"^" with 
the scalar product {A,B) ^ Tx{AB*). Write Wq = W^""* for the smallest *- 
closed subspace of IC'^'^^'^" which contains and P^"^ for the orthogonal 

projection of IK'^^X'^" onto W^"\ 

Theorem 1. Let C he a Wigner-Dyson or superconductor/ BdG sym- 
metry class. For u = 1, . . . ,r and n G N let A^'^''^^ G ^^SnxSn y^^-g^^ deter- 
ministic matrix such that for all v there exists a positive limit 

(2) lim -Tr[P^"^(A('^'"))(P^"VA('^'")))*l. 

n— >oo fi L u V / \ L, \ 

Then 

(i) As n — > CO, converges weakly to a random vari- 
able (i' = 1, . . . ,r). 

(ii) (Ti, . . . ,%■) has a centered joint Gaussian distribution, and there ex- 
ists a constant 7c > such that for all fijV = 1, . . . ,r 

(3) Cov(r^,T^) = Jim^^Re{Tr[P^"^(A(^'"))(P^")(^ 

Theorem 2. Let C be a chiral symmetry class. For u = l,...,r and 

n G N let A^'^'^^ G ]K"^" be a fixed deterministic matrix such that for all 
V there exists a positive limit (2). Assume furthermore that the sequences 
{p{n))n, {q{n))n [n =p{n) +q{n)] have the property that for all u there exists 
the limit 

*C,-:=Jl?i, Tr(P^ '{A^^ >)Ip(n),qin)), 
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ifC= BDI and AIII, and 

n— >oo fi ^ "^^ ' 

ifC = CII. Then 

(i) Asn^co, converges weakly to a random 
variable Ty (u = I, . . . ,r). 

(ii) (7i, . . . ,%■) has a Gaussian distribution with mean vector {tc,i, ■ ■ ■ ,tc,r) 
and covariance matrix given by (3) for a suitable 7c > 0. 

Remark 2.1. 

(i) 7c = i for C = A, All, 

7c = 1 for C = AI, AIII, B/D, Dili, C, 
7C = 2 for C = BDI, 
7C = 4 for C = CI, CII. 

(ii) For the B/D and C classes, Tr(^('''")yj"^) is in fact real. It follows 
immediately from the explicit description of the spaces Wc in Section 4 that 
the tc,v are real. 

We close this section with some technical notation that will be important 
for all subsequent sections. For sets SI, fi' write J-{^, fi') :={(/?: O —> O'}. If 
J) = {1, 2, . . . , A:}, S)' = {1, 2, . . . , n}, write T{k, n) := J^{n, 0'). If is even, 
write A4{i}) for the set of all pair partitions of Q, that is, partitions of 
all of whose blocks consist of exactly two elements. If SI = {1, 2, . . . , 2A;}, 
write M{2k) := M{n). For m G M{n) write JP"(m,0') for the set of all 
ip G S^') which are constant on the blocks of m. We sometimes write 
functions in the right operator fashion, that is, j<p in the place of (p{j). 

3. Extended Wick calculus. We will prove Theorems 1 and 2 via conver- 
gence of cumulants, hence of moments. To make this strategy operational, 
we need tools for the asymptotic evaluation of integrals of polynomial func- 
tions over the classical compact groups. These are available from recent 
work of Collins and Sniady [4]. To give a self-contained statement of these 
integral formulas, and to insert a sign factor into the symplectic integration 
formula which is missing in [4], we have to explain how these formulas are 
rooted in the invariant theory of the classical groups. As an extra benefit, 
this makes it possible to shed some light on how our tools relate to the 
Diaconis-Shahshahani integral. Let us start with this point, and focus on 
the conceptually simplest case, the orthogonal one. Consider a random vari- 
able Tn with values in 0„, distributed according to Haar measure. In order 
to prove via the method of moments that the random vector 

(Tr(r„),Tr(r2 ),..., Tr(r;)) 
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converges to a Gaussian limit as n — > oo, one has to evaluate matrix integrals 
of the form 

(4) / iTr{g)r {Trig')r ■ ■ ■ {Trig^r XoM 

[r G N, a = (ai, . . . , ar) E Ng]. As it turns out, the integral (4) has an inter- 
pretation in terms of the invariant theory of the group On- Specifically, let 
V be an n-dimensional complex vector space with scalar product (•,•) and 
symmetric bilinear form (3, such that 0(n,C) is the group of isometrics of 
P and 0„ = 0(n, C) n U„,, where U„ is the group of unitary matrices w.r.t. 
(•, •). Write p for the natural representation of 0(n,C) (or 0„) on V. On the 
tensor power V'^^ consider the representation = p®^ , which is given by 

/ fc \ k 
[^yApk{9)--=^Vjp{g). 
\j=i J j=i 

Define a representation of the symmetric group on V®^ via 

k 

\j=l J j=l 

Then the integral (4) equals 

(5) Tr(fjfc(s)|[y®fe]o„), 

where [1/®^]°" = {v ^ V®^ : vpk{g) = u Vg G 0„} is the space of invariants in 
Y^k representation p^, and s G has cycle type (l"i2"2 • • •r'^''), 

whence k = Y7j=ij(^j- See [16] for details. 

Now we turn to the problem whose solution in [4] is the basic tool for our 
present purposes. Given if,ip£ J^{k,n) (see notation at the end of Section 
2), what can one say about the integral 

(6) J t[9j<fi,H'^Onidg). 

If ej (j = 1, . . . , n) is an orthonormal basis of y, 99 G ^{k, n), write : = 
0^=1 Gjip- Then {e<^ : 99 G J-{k,n)} is a basis of y®'^, orthonormal w.r.t. the 
extension of f3 given by 

(fc fc \ fc 
j=l j=l I 3=\ 

Observe that (6) equals 
fc 

y Y[Piejipg,ej^(,)dg = J P{e^pk{g),e^)dg = f3(^J e^pk{g) dg,e^ . 



8 



B. COLLINS AND M. STOLZ 



Now, I— > / vpk{g)dg is the orthogonal projection of V®^ onto 
hence the integral (6) is the e^-coordinate of the projection of e^p onto the 
space of invariants. This clarifies in algebraic terms the family resemblance 
between the matrix integrals (4) and (6). We will pursue this idea a bit 
further and show that all ingredients of the orthogonal integration formula 
from [4] can be given transparent definitions in terms of the function spaces 
T{m,n). The link to the approach in [4], via Brauer algebras, is provided 
by the fact that if the underlying vector space has sufficiently large dimen- 
sion, the J-{m,n) describe a basis of the space of tensor invariants of the 
orthogonal group. Setting mo = {{1,2},. . .,{21- 1,2/}} e A^(2/), 

ei= 

</'GJP(mo,n) 

and 

ef2'={0is:sGS2z}, 

Weyl's First Fundamental Theorem of the invariant theory of the orthogonal 
group says that for k odd, [y^'^'jO" is trivial, whereas for k = 21 even, it is 
spanned by the S2;-orbit 9^ ^' . In fact, for n > I, this orbit is a basis of the 
space of invariants, and it is parametrized by M.{21) ([16], Corollary 3.11, 
Lemma 3.12). Specifically, for m G Al(2/), there exists Sra S S2; such that 

9lSm = X! ^'P- 

Set 

c(m,n) := f3{9lSra,0lSn). 

For n > I the matrix C := (c(m, n))n, ng_v((2«) is invertible, and its inverse 
(8) (Wgo„(m,n))n,,,:=C-i 

coincides with the orthogonal Weingarten function as defined in [4] (where 
the reader will also find a definition for the case n < I). To see this, note 
that if T is the 0„-equivariant isomorphism from End(y®') given 

by u{T{vi (g) ■ ■ ■®V2i)) = l3{u,V2®Vi®- ■ ■®V2i)vi®v^® ■ ■■®V2i-i (see [12], 
Section 4.3.2), then /5_b(p), as defined in [4], equation (15), equals {6iSm)T, 
where the Brauer diagram p equals the pair partition m if the upper and 
lower rows of the diagram are identified with even and odd numbers, respec- 
tively. 

Now, / eipPk{g)dg being an invariant, we see that it vanishes if k is odd, 
and that for k = 2l it can be written as J2m&M{k) 

Oimi9iSxn) ■ Then we have 
am= /5(e¥:.,6'/Sn)Wgo„(m,n), 

neM(2l) 
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meM{2l) 

/3(e<^, 6liS„)Wgo„ (m, n)l3{eiSm, e^) 

m,n 



Proposition 3.1. Let cp, ip e J^{k,n). Then 




In particular, the integral vanishes if k is odd. 

Turning to the symplectic case, let y be a 2n-dimensional complex vec- 
tor space with a skew-symmetric bilinear form a. Then there exists a ba- 
sis Cj {i = l,...,2n) of y such that a(ej,e„+j) = l [hence a(e„+i, Cj) = — 1] 
for i = 1, . . . ,n and a{ei,ej) = if |z — j\ ^ n. Introducing coordinates for 
g G Sp2„ via etg = Yd^igiiei, we obtain a{ei,ejg) = gi+n,j, a{ei,ej+ng) = 
gi+n,j+n, a{ei+n, ejg) = -gtj, a{ei+n, ej+ng) = -gi,j+n- Given functions (p, 
tp £ J-{k, n), a, P £ J-{k, {0, 1}) we can write for j = 1, . . . , /c 

9j{ip+na),j{ip+nl3) = '^{^j{^p+n{l-a))^ ^j{il>+nl3)9) i 

hence 

(9) =(-l)^^Sj^°y" a{e^+ri{i-a),e^+nf3Pk{9))dg 

= (-l)''~^^^"«(^e^+n(i-Q),y" e^+nf3Pki9)dg^, 

where the extension of a to the tensor power is defined as in (7). Since there 
exist no nontrivial invariants in odd tensor powers of V, we may assume 
that A: = 2/ is even. 

To describe a basis of [y^^'j^Pan^ ^e have to proceed slightly more care- 
fully than in the orthogonal case. An ordered pair partition of {1, ... , 21} is a 



hence 

P{ I e^Pki9)dg,e^ 



So we obtain the following 
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set m° = {{my,ni,) : = 1, . . . , /} of ordered pairs such that m := {{my,ny] : v - 
1,...,/} g7W(20. Set 

21 

where 

(10) ^(^'^'j)-|e,(,+„(i_,))(-l)-, if j = n,. 

For n>l, it follows from Weyl's First Fundamental Theorem and Section 
3.4 of [16] that 

{6mo:m° €M{21)} 

is a basis for [y^^ijSpa^^ where A^(2/) is identified with the set of all ordered 
pair partitions {{my^riy) : = 1, . . . , / } such that m-j^ < rij^ for all i^. 

Now let LP £ T{21, n) and a £ T{21, {0, 1}). We say that if + na£ T{m°,2n) 
if G !F{m,n) and {mya^nyo} = {0,1} for all u = 1, . . . ,1. Obviously, for 
a{6m'> , Sifi+na) not to vanish it is necessary and sufficient that + na £ 
!F{m° ,2n). In this case, there exist unique r],e such that 

a^(^u(r/,e, j),e<^+„Q^ /O. 

Specifically, i] = p and ue = 1 if, and only if, m^a = 0. For these rj, e we have 

(11) a{9mo,e^+na) = a{^v{'q,£,j),e^+no^ =±1. 

To determine the correct sign, let us disregard the sign in (10) for a moment. 
Any = 1, . . . , / yields a factor of the form a(ej+n, ej)a(ei, Cj+n) = —1 in (11). 
So we have a factor (—1)' regardless of the specific choice for e. Now, we get 
an extra minus sign from (10) if ue = 1, hence rriya = 0. So we see that 

(12) (ll) = (_iy+#{i':'n."=0}_ 

It remains to compute (9). To this end, observe that k = 21 being even, 
a is a symmetric bilinear form on V®^ . So the corresponding steps of the 
orthogonal case are available. Applying (12) to <^ + n(l — a) and ip + n(i and 
keeping in mind that k = 21 is even, we see that (9) equals 

(13) {-lfir-i"=^} 

X J2 Wgsp2jm°,n°)(-l)#^^^™''"=^>+#^'^^"'^^=°> 

m°,n''eA^{2/) 

X l.^(m° ,2n){'P + na) l^(no _2n) ( + "/?) , 
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where we have genericahy written 

(14) m° ={{m^,n^):u = l,...,l}, n° = {{ml,n'^) = 1, . . . ,1}. 

Observe that all summands vanish unless the prefactor (13) equals (— 1)'. 
So the overall sign factor is 

Summing up, we have the following: 

Proposition 3.2. Let ip, £ J^{k,n), a, P £ T{k, {0,1}). Then 

k 

= Wgsp^Jm°,n°)(-l)#{'^^'"^"=°>+#^^^'<^=°^ 

m°,n°eA1(fc) 

X l.^(m» ,2n) i'P + l.^Cn" ,2n) (ip+nP), 

where we adopt the notation of (I4) above and identify A4{k) (k = 2l) with 
the set of all ordered Ui,) : = 1, . . . , /} such that ruy < 

for all v. 

Remark 3.3. Since gj = Jg for g G Sp2„, Proposition 3.2 covers also 
the case of complex conjugates in the integrand. 

Finally, let us record the unitary integration formula. 

Proposition 3.4. Let ip,ip,ip',ip' € J^{k,n). Then 

/ n ^JVJ'/'^jV' jV'' ^^nidu)= '5('/^'Cr<^')'^(^>7-^')Wgu„(T"cr"^), 

where the unitary Weingarten function Wgy^ is the inverse for convolution 
on the symmetric group S^ of the function a ^ 72#cycles{fT) _ 

In what follows, there will be no need to work with the definitions of the 
various Weingarten functions. It suffices to know their fundamental asymp- 
totic properties: 

Proposition 3.5. 

/icN inr / \ ( n~''(l + 0(n~^)), ifa = idk, 

^^U"(") = |o(n-(^+i)), otherw^^. 
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For k = 2l, we have 



(17) Wgsp^Jm;,m5) 



n 



(l + 0(n-i)), if mi=m2e M{k), 



0{n '•'^^^), otherwise. 



Proof. [4], Theorem 3.13. □ 

in) 

4. The proofs. Recall that we write for the *-closed K-linear hull 

of S(f^ in IfC'^"x<5"j and for the orthogonal projection of IC^"^^'^" onto 

W^""* with respect to the scalar product {A^B) i-^Tr(A_B*). Then obviously 
Re(Tr(^('^'")yj"^)) = Re(Tr(P(^"^(A('''"))yc^"^)), so we can restrict our atten- 
tion to taken from W^"\ Observe that the map W^"^ M,yl(^'") ^ 

Re(Tr(^('''"Vc^"^)) is M-linear. By Cramer-Wold (see, e.g., [2], Theorem 
29.4), this implies that it suffices to prove convergence to a Gaussian limit of 
variance lim„ ^ Tr(A('^'")(A('^'"))*) = lim„ 2£ Tr[P^"^(^(^'"))(P^"^(A(^'")))*] 
for an individual v. Recall that for a Gaussian distribution, the third and 
higher cumulants vanish. Denote the rth cumulant of Re(Tr(^('^'"') V^^"^)) by 
C^''\ It follows from [4], Theorem 3.16, that for r > 3, d"'' = o(C^"^), as 

n ^ oo. So it suffices to show that Cg^^ that is, the variance, converges to 
the correct finite limit. We verify this in a case- by-case fashion. 

4.1. Class A. We have $^^(5) = g{g G U„) and W^^'^ = C"^". By abuse 
of language, we also write g = g^""^ for a random element of Un, with distribu- 
tion given by Haar measure. Haar measure on U.„ being invariant under mul- 
tiplication by scalars z G C with |z| = 1, we see that E(ReTr(A(")5("))) =0 
for all G W^^ (n G N). Now, 

E((ReTr(A(")5(")))2) =Ef ( ^ ^ ; ^ g ) 

(18) J 

= -E(Tr(A('^)5('^))TV(^H(7W)). 



The second equality is evident from Proposition 3.4. Now, (18) equals 

n . n 

\ / 9jiW^ ^vSdg) = ^WguJidi) 



(n) (n) 

by Proposition 3.4. By Proposition 3.5, this is asymptotically equivalent to 

— Tr(A(")(A("))*). 
2n 
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4.2. Class AI. To simplify notation, we will from now on suppress de- 
pendence on n in most cases. We have ^Ai{g) = 99' id G U^), Wai = {A G 
C"^" : A' = A}. Recall from Section 2 that we write V = = F^j^ for a 

(n) 

random element of , distributed according to the push-forward of Haar 

(n) 

measure under . Continuing to abuse language, and to write g for a Haar 
distributed element of Un, we have V = gg' . As in the A case, one sees that 
E(ReTr(^y)) =0 and 

E{{ReTT{AV)f) = ^E{Tt{AV)Tt{AV)). 

Now this equals 



i,j,k,i' ,j' ,k 

which is asymptotically equivalent to 



h X] o.jiO'j'i' I 9ik9jkgi'k' 9j'k' dg, 



hence to 

-Tr(A(")(A('^))*), 
n 

since A' = A. 

4.3. Class AIL Here $aii(5) = ^^V J (ff G Uan) and VFaii = {A G C2"x2n . 
{Ajy = -AJ}. As in the A and AI cases, E(ReTr(^(")y("))) = and 

E((ReTr(^("V(")))^) = iE(TV(^("V(''))Tr(AWy W)). 

Write J = (Jij). One has TiiAV) = TT{AgJ'g'J) = -Tr{BgJg') with B := 
J^. Now, 

i 1 TV(i?5Jff')Tr(SgJ<7')c?5 

= ^ X! bijJkih'j'Jk'V j 9jk9ii9j'k' 9i'i' dg. 

i,j,k,l,i',j' ,k' ,1' 

This is asymptotically equivalent to 



iWgu2„(id2)( Y bijJkibijJki + bijJkibjiJik 

\i,j,k,l i,j,k,l / 



hence to 



^ (Tr(SS*) Tr( JJ') +Tr(SS)Tr(j2)) = -^—Tr{AA*)Tt{JJ'). 



2(2n)2' ' / V / ' V ; V ^2n)2 
So we have obtained that E((ReTr(A(")y(")))2) ~ ^ Tr(A(")(yl("))*). 



14 B. COLLINS AND M. STOLZ 

4.4. Class AIII. In this case one has ^Aiii{g) = glp,qg* Ip,qi 9 G U„ (n = 
p + q) and VFaiii = {A G C"^": A/p,, = (A/p,^)*}. For ^ = ^(") G W^l"^ 
we have Tr(AF(")) = Tr(^c//p,5£f*/p,g) = TriBglp^gg*) with B := Ob- 
serve that this is a real random variable. Writing Ip q = (lij), we have that 
E{TT{BgIp^gg*)) =J2i,j,k,i^ijhi J gjkgudg is asymptotically equivalent to 



l;J2buhk = l^TT{B)TT{Ip 

(19) 



n ^ n 



p{n)-q{n) p{n) - q{n) 
TV(5) = Tr(^^ 'Ip{nU{n))- 



n n 

On the other hand, 



E(Tr(Sg/p,q5(*))2 = kjbi'j'Ikih'i' 9jkgj'k'gu9i'i' dg 



i,j,k,l,i' ,j' ,k' ,1' 

is asymptotically equivalent to 



Wgu„(id2)( kih'i'hkh'k' + Y kvh'ilvihvy. 



hence to 

^^((Tri?)2(Tr/p,,)2+Tr(i?2)TV(/^_^)). 



Now (19), the assumptions of Theorem 2, and the definitions of B and Waiii 
imply that as n — > oo, 

V(Tr(S(")5(")/p,g(5("))*)) = V(Tr(A("V("))) ~ i TV(A(")(^("))*). 

4.5. Class BID. Here «>b/d(5) = 5,5 e 0„, and Wb/d = IR""""- By in- 
variance of Haar measure for A G Wg/D o'^^ has E(Tr(j45()) = 0. On the 
other hand, by Proposition 3.1 

E(TV(yl<7)') = E "^J-^'^^ / 5,i5«fcd5 = WgoJ{{l,2}},{{l,2}})E4-, 

i,j,k,l i,j 

which is by Proposition 3.5 asymptotically equivalent to ^ Tr{A^'^\A^'^^y). 

4.6. Class BDI. In this case one has <^bdi(5) = glp,qg'lp,q, g G 0„, 
n = p + g and Wrdi = {A G M"^" : (A/p,^)' = A/p,J. Then Tr(^y(")) = 
Tr{AgIp^gg'Ip^g) = Tr:{BgIp^gg') with B := Ip^qA. Writing again Ip^q = (lij), 
one obtains 

E{Tr{BgIp^gg')) = ^ kjhi J gjkgudg, 

i,j,k,l 
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which is asymptoticahy equivalent to 

(20) - ^ biJkk = - Tr(i?) Tr(/p,,) = ^^""^ ~ ^^""^ T¥(A(")/p(„),,(„)). 

i,k 

On the other hand, 

E((Tr(Sg/p^gg'))^) = X! bijbi'jJkih'i' J 9jk9u9j'k'9i'i' dg. 

i,j,k,l,i' ,j' ,k' ,1' 

By Propositions 3.1 and 3.5, this is asymptoticahy equivalent to 



(21) 

n 



li X! biibin/Ikkh'k' + X! bii/bifilkk'h'k + X! bijbijikilkij 

\i,i',k,k' i,i',k,k' / 

= -^2((Tri?)2(Tr/,,j2 +Tr(S2)Tr(/2j +Tr(Si3')T^-(W;,,)) 
= ^((Tr5)2(p - qf + 2nTV(55')), 

where the three sums in (21) correspond to the pair partitions {{1, 2}, {3, 4}}, 
{{1,4}, {2, 3}} and {{1, 3}, {2, 4}} of {1, 2, 3, 4}. Comparing with (20), we see 
that V(Tr(AWy("))) ~ f T¥(^(")(A("))'). 

4.7. Class Din. Here one has <&diii(s') = 9J'9'J,9 G 02n, and Wdiii = 
{A G M2«x2n : (^j)' = - A J}. Then Tr(^Wy(")) = Tr(ylc/JVJ) = Ti{BgJ'g') 
with S := Jyl. 

'E{TT{Bgj'g')) = ^ 6ijJ/fc y gjkgudg = ^"^buJkk = 

by definition of J. On the other hand, 

E((Tr(55rJV))^) = ^ bijJikbiijiJvk' j gjkgii9j>k'9i>i'dg. 

i,j,k,l,i' ,j' ,k' ,1' 

As in the BDI case, one is led to sums corresponding to the pair partitions 
of {1,2,3,4}. The sum corresponding to {{1, 2}, {3, 4}} vanishes, because it 
involves diagonal terms of J. The sum corresponding to {{1, 3}, {2, 4}} is 
(2^ E^,J,k,l hjkjJikJik = Tt{BB') TV( JJ') = ^ Tv{AA'). The last sum 
^^Ji^E^,^',^'bu'Jwb,,Jn = ^Tr(i32)Tr(j2) = (-Tt{AA')){-^). Sum- 
ming up, 

E((Tr(A("V(")))^)~-T¥(A(")(^("))'). 
n 
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4.8. Class C. Here we have ^c{g) =9,9^ Sp2„, and Wq = H"^". Since 
jjnxn -g giQggd under multiplication, we see that Tr{Ag) is real for 



(22) A=( ^ 



wnxn 



and the same argument applies to classes CI and CII below. Invariance of 
Haar measure yields K{Ti{Ag)) = 0. Furthermore, from Proposition 3.2 we 
obtain 

E((Tr(A("V(")))^) 
= E((TY(AW5))') 

= ^ ^ 0,j+nt,i+ns(lj'+nt' ,i'+ns' J 9i+ns,j+nt9i'+ns' ,j'+nt' dg 
s,t,s',t' 
e{l,...,n} e{0,l} 

= WgSp,,J{(l, 2)}, {(1, 2)}) ^(-1)''+* aj+nt,i+nsaj+n{l-t),i+n{l-s) 



^ji^j+riyi+n O-j+nA+nO-ji CLj+n,iOjj,i+n CLj, 1+11(1 j+n,i 



n ^ 

-(Tr(X(X)') - Tr(y(-F)')) = - Tr(XX* + YY* 

n n 



-TrM(")M(")) 
n 



4.9. Class CI. We have $01(5) =S'^,nS'~^^,n(5'eSp2„) and Wci = {^G 
"X":(A/„,„)*=^/„,„}. Set 

In 
In 



K . — JIn,n — In,nJ " 



and write K = (nij). Since gJ = Jg for g G Sp2„, we have for any A € Wc\ 
that 

= Tr(A5r/„,„5r~^/„,„,) 

= TV(A(7/„,„(-J)(7'J/n,n) 

= Tr(A5i^5'K)=Tr(BgAV) 

i ,i,A;,Ze{l,...,n} g,r,s,te{0,l} 

where iiTA = : B = {hij). It follows easily from Proposition 3.2 and the sym- 
metry of K that E(Tr(SgKc?')) = 0. Turning to ¥.\{TT{BgKg')f), we apply 
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Proposition 3.2 to 

^ ^ ^ ^ bi-\-nq,j+nrbi'-\-nq' ,j'+nr' ^k+ns,l+nt^k'+ns' ,l'+nt' 
i,j,k,l q,r,s,t 
i',j',k',l' q',r',s',t' 



The ordered pair partition {(1,2), (3,4)} leads to summands of the form 

~lj Oi+ra(l-r),i+nr''i'+n(l-r'),i'+nr''^^fc+ns,fc+7i(l-s)'*fc'+ns',fe'+7i(l-s')' 

being symmetric, the summands come in it pairs, and the total sum 
vanishes. To {(1,3), (2,4)} correspond summands of the form 

{~'^)'''^^~^^~^^bi+nqj+nrbi^n{l-q)j+n{l-r)''^k+ns,l+nt>^k+n{l-s),l+n{l--t]- 

In view of the specific form of K, for such a summand not to vanish it is 
necessary that k = I and s ^t. Hence the contribution from {(1,3), (2,4)} is 

2?^— 2 ^ {~^y~^^~^^biJ^nq,j+nrbi-\.n{l-q),j+n{l-r) 

(23) ^ 

Note that we have obtained the factor 2n since we have suppressed the sums 
over k and s. Writing A as in (22) above, we have 

^={x y)' 

and since AIn n was assumed to be hermitian, we have X = X and F = y'. 
This means bi+nj+n = bj+n,i+n = -hj and 6jj+„ = 6j+„,j- = 6j,j+n. Hence, 
(23) equals 

-(Tr(XX*) +Tr(yy*)) = -Tr(ylyl*). 

n n 

Finally, the ordered pair partition {(1, 4), (2, 3)} yields summands of the 
form 

(~l)''^'^^'^^*^i+ni},j+nT^j+n(l-r),i+n(l-(j)'^fc+nsJ+ntKi+n(l-f),fe+n(l-s)5 

and an analogous argument shows that the contribution of {(1,4), (2,3)} is 
2 2 

^ ^ bi^j-\-nbj^i-\-n ~\~ biJ^n,jbj-\-n,i bi^jbjj^n,i+n bi-^n,j+nbj^i — Tr(j4j4 ). 
ij 

Summing up, E((Tr(A(")y W))^) = E((Tr(B5i^5'))^) ~ | Tr(A('^)(^("))*). 
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4.10. Class CII. We have ^cii{g) = 9'^p,qg~^'^p,q {g G Sp2„, n = p + 
Wen = G EI"X" : (^p,,)* = Alp^q}- Write 

K —IT - ( ^ ~h,q\ _r / 



and Kp^g = (nij). Fix 



Then 

Tr{AV) = Tv{Aglp,gg~%^g) = TviAgIpJ-J)g'JIp^g) = -Tt{BgKp,gg' 
where (bij) = B := Kp^qA. Note that 

(24) B=( ^P''^^ ~^P^i-^ j _ 

Since Alp q is assumed to be Hermitian, we obtain the useful relations 

(25) Xlp^q = (Xlp^q)* = Ip^gX* , YIp q = —{YIp,q)' = —Ip qY'. 

Now, ¥.{-TT{BgKp^qg')) equals 

^ ^ ^ ^ bi+nq,j+nrHk+ns,l+nt j 9j+nr,k+nsgi+nq,l+ntdg. 

hj,k,l€{l,...,n} q,r,s,t£{0,l} 

The leading order contribution comes from the case j = i, k = 1, q = {1 — 
t = (1 — s), hence is asymptotically equivalent to 

(26) ^ ^( — l)''''''^6i-|_n(l-r),i+nr^fc+ns,fc+n(l-s) 



1 ' " 



= -- \^Y.>^w^ Civ(/p,g(x + X))) 

= 2 (^(")7(-)) Tr(AHT,(.,,,.)). 
Turning to the second moment, we have that "^{{^{BgKp^qg'))^) equals 



^ ^ ^ ^ bi+nq,j+nrbi'^nq',j'+nr'^k+ns,l+ntf^k'+ns',l'+nt' 
i,j,k,l q,r,s,t 
i',j',k',l' q',r',s',t' 



^ / gj+nr,k+nsgi+nq,l+ntgj'+nr',k'+ns'gi'+nq',l'+nt' dg. 
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First of all, the leading order contribution corresponding to the ordered pair 
partition {(1, 2), (3, 4)} is asymptotically equivalent to 

~ \ ^) "i+n(l-r),j+nr'^j'+n(l-r'),i'+nr' 

i,k,i' ,k' r,s,r',s' 

^ l^k+ns ,k+n(l-s)'^k' +ns' ,k' +n{l-s') J 

which coincides with the square of (26). So, asymptotically, the variance will 
come from the remaining ordered pair partitions. Next, from {(1,3), (2,4)}, 
we get summands of the form 

{~'^y'~^^~^'^~^^bi+nqJ+nrbi-^.n{l-q)J+n{l-r)'^k+ns,l+ntl^k+n{l-s),l+n{l-t)■ 

For such a summand not to vanish, it is necessary that k = I and s ^ t. 
Note that in this case the product of the k terms is —1. So we obtain a 
contribution 

'^^~2 ^y^.^y^X~^y^^^i+nq,j+nrbi-\-n(l-q),i+n(l-r) 

— ^ ^ bijbn-\-i,n+j ^i,n+j^n+i,j 
i,j 

= ^(T¥((/p,gF)(/p,,y)0 +Tr((/p,,X)(/p,,X)')) 

= -(Ti(XX*) + Tr(YY*)) = - Tr( AA*). 
n n 

Finally, as to {(1,4), (2,3)} one argues similarly to obtain a contribution 

+ Tr((/p,,X)2)+Tr((/p,,X)2)), 
and one invokes (25) to see that this in fact equals ^ Tr(A(") (A^"))*). 
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